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Within a minimal model, we present analytical expressions for the eigenstates and eigenvalues 
of carriers confined in quantum rings in monolayer and bilayer graphene. The calculations were 
performed in the context of the continuum model, by solving the Dirac equation for a zero width 
ring geometry, i.e. by freezing out the carrier radial motion. We include the effect of an external 
magnetic field and show the appearance of Aharonov-Bohm oscillations and of a non-zero gap in 
the spectrum. Our minimal model gives insight in the energy spectrum of graphene-based quantum 
rings and models different aspects of finite width rings. 
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I. INTRODUCTION 

The investigation of low-dimensional solid state devices 
has allowed the direct observation of quantum behav- 
ior in electron systems. These effects arise due to the 
confinement of carriers in structures that constrain their 
movement along one or more directions, such as quantum 
wells, quantum wires and quantum dots. One important 
class of such low-dimensional systems are quantum rings, 
in which a particular type of confinement together with 
phase coherence of the electron wavefunction allows the 
observation of effects such as the Aharonov-Bohm^ and 
Aharonov-Casher~ effects. Quantum rings have been ex- 
tensively studied in semiconductor systems, both experi- 
mentally and theoretically^ and are expected to find ap- 
plication in microelectronics, as well as in future quantum 
information devices. 

In this paper we present analytical results for the eigen- 
states and energy levels of ideal quantum rings created 
with graphene and bilayer graphene. Graphene is an 
atomic layer of crystal carbon which has been the tar- 
get of intense scrutiny since it has been experimentally 
produced^—. Part of this interest stems from the unusual 
properties of carriers in graphene, caused by the gap- 
less and approximately linear carrier spectrum, together 
with possible technological applications, such as transis- 
tors, gas sensors and transparent conducting materials 
in e.g. photovoltaics. Additionally, it has been found 
that two coupled graphene sheets, also known as bilayer 
graphene displays properties that are distinct from sin- 
gle layer graphene as well as from graphite. The car- 
rier spectrum of bilayer graphene is gapless and approx- 
imately parabolic at the vicinity of two points in the 
Brillouin zone^. In particular, the spectrum is strongly 
influenced by an external electric field perpendicular to 
the bilayer, leading to the appearance of a gapi .. The 
high quality of the single layer and bilayer graphene sam- 
ples that have been obtained, together with the large 
mean free path of carriers suggests that phase coher- 
ence effects may be observable in these systems. Re- 



cently, graphene-based quantum rings produced by litho- 
graphic techniques have been investigated on single-layer 
graphen o 11 ' 12 . These systems have been studied theoret- 
ically by means of a tight-binding model, which does not 
provide straightforward analytical solutions for the eigen- 
states and eigenvalues^—. For bilayer graphene also it 
was recently showr>i& that is possible to electrostatically 
confine quantum rings with a finite width. The energy 
spectrum was obtained by solving the Dirac equation nu- 
merically. 

In this paper we present a toy model that allows for 
analytical expressions for the energy levels of quantum 
rings in single layer and bilayer graphene. This model 
permits the description of several aspects of the physics 
of graphene quantum rings without the additional com- 
plications of edge effects and finite width of the quantum 
ring. We are able to obtain analytical expressions for 
the energy spectrum and the corresponding wavefunc- 
tions, the persistent current, and the orbital momentum 
as function of ring radius, total momentum and magnetic 
field, which can be related to the numerical results ob- 
tained by other methods. 

The paper is organized as follows: in section II we 
present the theoretical model and numerical results for 
quantum rings in single layer graphene. Similar results 
for bilayer graphene are given in Sect. III. Section IV 
contains a summary of the main results and the conclu- 
sions. 



II. SINGLE LAYER GRAPHENE 
A. Model 

The dynamics of carriers in the honeycomb lattice of 
covalent-bond carbon atoms of single layer graphene can 
be described by the Dirac Hamiltonian (valid for E < 0.8 
eV). In the presence of a uniform magnetic field B = 
Boe z perpendicular to the plane and finite mass term 
A, which might be caused by an interaction with the 
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underlaying substrat o 10 ' 19 The Hamiltonian in the valley 
isotropic form is given by-i^: 



H = v F (p + eA).er + rAa z 



(1) 



where r = +1 corresponds to the K point and r = — 1 
to the K' point, p is the in-plane momentum operator, 
A is the vector potential and v F w 1.0 x 10 6 m/s is the 
Fermi velocity, and cr = (u x ,a y ,a z ) is the pseudospin 
operator with components given by Pauli matrices. The 
eigenstates of Eq. ([1]) are two-component spinors which, 
in polar coordinates is given by 



<t> A {p)e im * 
i(j) B {p)e i{m+1) 



(2) 



where m is the angular momentum label. We follow 
the earlier very successful approac h 20 ' 21 of ideal one- 
dimensional (ID) quantum rings in semiconductors with 
spin-orbit interaction where the Schrodinger equation 
was simplified by discarding the radial variation of the 
electron wave function. Thus, in the case of an ideal ring 
with radius R, the momentum of the carriers in the ra- 
dial direction is zero. We treat the radial parts of the 
spinors in Eq. ([2]) as a constant 



cj) A (R)e lm ' t ' 
i(j)B{R)e l{m+1)4 ' 



(3) 



Because the radial motion is frozen in our model 
there will be no radial current and the persistent cur- 
rent will be purely in the angular direction. By solving 
H'^{R 1 4>) = E^(R 1 <j)) and using the symmetric gauge 
A = (0, Bqp/2, 0), we obtain 



(m + 1 + P)<t>B(R) = (e - tS)^a(R), 
(m + /3)MR) = (e + T5)ct> B (R), 



(4) 



where the energy and mass terms are written in dimen- 
sionless units as e = E/Eq, S = A/Eq with Eq = hv F /R. 
The parameter (3 = (eBo/2h)R 2 can be expressed as 
(3 = $/$o with $ = ttR 2 Bo being the magnetic flux 
threading the ring and $o = h/e the quantum of mag- 
netic flux. The homogeneous set of equations ((3]) is solv- 
able for the energies 



e = ±V(™ + /3 + l)(m + /3)+(5 2 . 
This energy can also be written as 



= ±-\/ (in — m_)(m — m_|_) 



where 



m± 



-(/?■ 



± 



5 2 . 



(5) 



(6) 



(7) 



Note that the energy spectrum for an ideal single layer 
quantum ring for both K and K' points are the same. 
For \5\ > 1/2 we have that m + = m*_ is complex and 



e is real for any value of /3. In the region — | 



<5< I 



the energy is real, except for to_ < m < m + , when the 
energy is complex. For the gapless case, i.e. S ~ 0, we 
have m+ = —f3 and m_ = — (3 — 1 and the energy is 
real when m < — /3 — 1 or when m > —j3 and imaginary 
otherwise. 

The wavefunctions are eigenfunctions of the total an- 
gular momentum operator given by the sum of orbital 
angular momentum L z and a term describing the pseudo- 
spin S z 



Jz — L z + hS z 



(8) 



where S z — (l/2)er z , with a z being one of the Pauli 
matrices and the eigenvalues of J z operator become 
[m+(l/2)]ft. 

The current is obtained using jx,y = "f[^V Ii9 ^]. The 
total angular current j = vp^^ cr^tl'] can be calculated 
using the fact that = £(0)er. y , where 



e"^ 
e i4 



(9) 



The current for the electrons in the K- valley becomes 



]i< = v F ((p A (p B 



The total current is then given by j = jx 4 
radial part of the two spinor components are 



3K> 



MR) = i, <Mi?) = 



e + rS' 



(10) 
The 

(11) 



Notice that the radial current can be calculated using 
j r = vf[iI>' £,(4 > )< J x'ifi) which leads to the the following re- 
lation, 



j r = iv F {(j) A (j)B - <I>bM, 



(12) 



where, in the case of our ideal ring we have j r = 0. From 
Eqs. (TTUJ) and ([PI]) , one can find the following expression 
for the total angular current of a single layer quantum 
ring 



J = 



Av F e(m + /3) 



(13) 



e 2_ S 2 ■ 

One can rewrite Eq. (|13[) in the following form 

A = A +(— ) l 2 ( m + P)(c 2 + 62 )-(t 2 - S2 ) 
v F { dj3 >K ^dj3 >K ' e(e 2 -(5 2 ) 

(14) 

Since for the ground state energy y/S 2 — 1/4 < e < 6 
and — 1/2 < m + (3 < the last term in Eq. (13) is much 
smaller than the derivatives of the energy with respect 
to the flux ($) and oscillates around zero. Note that in 
ID semiconductor rings the current is exactly given by 
dE/d§, which is thus different from graphene where we 
have approximately 



J_ 

v F 



- de s 
•dj3' K 



, d p)K> 



(15) 



with (3 = ^/f&Q. 



3 



B. Results 



The energies as function of ring radius R are shown in 
Fig. [TJ for A = 50 meV, with — 10 < m < — 1 (magenta 
curves), 1 < to < 10 (blue curves), and m = (green 
curves). In the absence of an external magnetic field, 
the energy is given by E = ±^Jm(m + l)(hv p / R) 2 + A 2 
and the energy branches have a 1/R dependence and 
approach E — > ±A for very large radii. Note that for 
m = and m = — 1 the energy E = ±A is indepen- 
dent of R and all branches are two-fold degenerate. For 
non-zero magnetic field (B = 3T), the right panel shows 
that the branches have an approximately linear depen- 
dence on the ring radius for large R, in particular we 
have E ~ ±^{aR) 2 + A 2 , with a = VFeBo/2. For small 
radii, E ~ ±hyjm(m+ 1)/R and all branches diverge as 
1/R, except for m = and m = — 1. In those cases when 
ii^Owe have for to = the result E = ±\/A 2 + ahvp, 
while m = — 1 gives E = ±\/A 2 — ahvp. 

Figure [5] presents results for the energy as function 
of total angular momentum index to, for A = 50 meV, 
R = 50 nm and for three different values of magnetic 
field, namely B$ = — 5 T (diamonds), Bo = 0T (circles) 
and Bq = 5 T (triangles) . Notice that for a given Bq the 
electron energy obtains a minimum for a particular to, 
i.e. for B Q = (5 T, -5 T) it is m = (9, -10). In fact 
it is given by m = — (&/§o + 1/2) and is independent of 
A. 

The energy levels as function of the external magnetic 
field are shown in Fig. [31 for a quantum ring with (a) 
6 = 1/2, (b) 5 = 3/8, (c) 5 = 1/4 and (d) 6 = with 
R = 50 nm for — 10 < m < — 1 (red curves), 1 < to < 10 
(blue curves), and to = (green curves). The magnetic 
field dependence of the spectrum becomes evident if one 
rewrites Eq. (5) as e 2 - [(m + $/$ ) + 1/2] 2 = S 2 - 
1/4. Thus, for the special case of <5 = A/Eq = 1/2 the 
gap is zero and the energy levels are straight lines given 
by e = ±(m + l/2 + $/$o)- The energy spectra for 
5 > 1/2 resemble those found earlier by Recher et alJ£ 
in the case of a finite width graphene ring with infinite 
mass boundary conditions. An enlargement of Fig. [3] 
around E = is shown in Fig. 2) The spectrum has 
an interesting magnetic field dependence with decreasing 
5. For (5 = the double degeneracy is restored at E — 
0. This behavior can be easily illustrated by considering 
to = 0. The energy in this case is e = ±^ f3([3 + 1) + S 2 
which for 5 — 1/2 becomes e = ±(/3 + l/2) while for 5 = 
it is e = ±y/p(/3 + l) and thus e ~ ±V/J for ~ 0. 

In Fig. [S^a) the energy spectrum is plotted vs mag- 
netic field for (5 = 2. Where, the energy has a hyperbolic 
dependence on the applied magnetic field with minima at 
$/$o = -to - 1/2 and a gap of Ae = 2y/5 2 - 1/4. The 
exact location of the transitions (orange dots) and the 
location of the minima points (yellow dots) in the energy 
spectrum is clarified in Fig. [S](b). The dependence of the 
energy levels on the gap parameter A is shown in Fig. [51 
for zero magnetic field (left panels) and B^ = 1 T (right 



FIG. 1: (Color online) Energy levels with m = — 10, 10 of 
single layer graphene quantum ring as function of ring radius 
R for Bq — T (left panel) and Bo = 3 T (right panel) when 
the mass term is A = 50 meV. 



panels). When to > (upper panels) and to < (lower 
panels). For Bo = T the energy levels are two fold de- 
generate where E{m) = —E(—m — 1). When a magnetic 
field is applied an energy gap is opened (see right-bottom 
panel in Fig. 6). Notice also that the m = —2 level only 
exists for A > E /2, i. e. for A < Eq/2 there is no real 
energy solution when to = — 2. 

The corresponding ground state expectation values for 
the operators in Eq. (8) are plotted as function of the 
magnetic field in Fig. [TJb) for both K (black dashed 
curve) and K' valley (black dash-dotted curve). Notice 
that for the if -valley < L z >~ mh and < S z >~ h/2 
whereas in the if '-valley < L z >~ (m+l)H and < S z >~ 
—H/2. Thus for both the K valley and the K' valley 
< J z >~ [to + (1/2)]?! which is approximately quantized 
and on the average its value decreases linearly with the 
applied magnetic field. 

The angular current density for a single layer graphene 
quantum ring is shown in Fig. [5{c). Note that the con- 
tribution from the if -valley jjc (Fig. UK a)) an d the K'- 
valley (Fig. (HUb)) are not the same, they have oppo- 
site sign and oscillate in phase around a nonzero average 
value —rvp/A. The reason is that if for given energy we 
have electrons in the if -valley, the corresponding parti- 
cles in the if '-valley will behave as holes. The persistent 
current is a sawtooth shaped oscillating function of the 
magnetic field with period ABq = ^q/ttR 2 . This be- 
havior is quantitatively very similar to those found for 
the standard Aharanov-Bohm oscillations in metallic and 
semiconductor quantum rings. 
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FIG. 2: (Color online) Energy levels of a single layer graphene 
quantum ring as function of the quantum number m for Bo = 
-5, 0, 5 T with A = 50 meV and R = 50 nm. 



III. BILAYER GRAPHENE 
A. Model 

In the case of bilayer graphene the Hamiltonian in the 
vicinity of the K point, is given by^ 



H = 



( t\J\ 7T t 

Ti-1" tU\ 

t tU 2 tt 

\ TT tU 2 J 



\ 



t 



(16) 



where r = ±1 distinguishes the two K and K' val- 
leys. £ ~ 400 meV is the interlayer coupling term, 
7T = vp[(p x + eA x ) + i(p y + eAy)], U\ and U 2 are the 
potentials, respectively, at the two graphene layers. 
Here we do not include any mass term because the 
gate potential across the bilayer is able to open an 
energy gap in the spectrum^ The eigenstates of the 
Hamiltonian (jT5J) , are four-component spinors ty(r,(f>) = 
[Mp)e im *jJMp)e i(m ~ 1) * > 0c(p)e im *, i^ D {p)e l{ - m+1 ^] 
(see Ref. |22j)- Following our earlier approach for an 
ideal ring with radius R, the wave function becomes: 



<F = 



(l)c{R)e im 't> 
\ i<t> D (R)e^ m+1 ^ / 



(17) 



We use the symmetric gauge and obtain the following set 
of coupled algebraic equations 

-(e - tui)4>a{R) -{m + p- l)<p B (R) + t'4> c {R) = 0, 
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FIG. 3: (Color online) Electron and hole energy levels of a 
single layer graphene quantum ring as function of external 
magnetic field B for (a) 8 = 1/2, (b) 8 = 3/8, (c) 5 = 1/4 
and (d) 5 = with R — 50 nm, and total angular quantum 
number —10 < ra < —1 (red curves), 1 < m < 10 (blue 
curves) and m — (green curves). 



.'=1/2 (5 = 3/8 5=1/4 6 = 




FIG. 4: (Color online) The same as Fig. 3, for m = — 1 (red 
curves), m = 1 (blue curves) and m = (green curves) and 
different values of dimensionless mass term 8. 



t'MR) - (e - ru 2 )(j> c {R) + {m + /3 + l)<f> D (R) = 0, 
(m + p)<t>c(R)- (e-Tti 2 )Mii) = 0. (18) 



(m + P)^ A {R) + (e - T Ul )<f> B (R) = 0, 



where, t' = t/E and U\ 2 = Ui i2 /E are in dimension- 
less units. After some straightforward algebra we obtain 
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FIG. 5: (Color online) (a) Electron and hole energy levels of 
a single layer graphene quantum ring as function of external 
magnetic field Bo for S — 2 and R = 50 nm. (b) An Enlarge- 
ment of the region which is shown in (a) by a rectangle. 
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the following polynomial equation that determines the 
energy spectrum 

(e - t Ui ) 2 [(e - tu 2 ) 2 - (m + /3)(m + + 1)] 

-(m + 0)(m + p- l)[(e-Tu 2 ) 2 - (m + /3)(m + p + 1)] 

-(e-rui)(e-Tit 2 )t /2 = 0. (19) 

After introducing the average potential w = (u\ +U2)/2 
and half the potential difference 8 = (u\ — ui)j2 we can 
rewrite this quartic algebraic equation in a more compre- 
hensive form: 

s 4 -2s 2 [(to + /3) 2 +8 2 + (t') 2 /2] 
+4:ST8(m + P) + (to + /3) 2 [(m + /3) 2 - 1] 



FIG. 6: (Color online) Lowest energy levels of a single layer 
graphene quantum ring as function of the mass term A with 
Bo = T (left panels) and Bo = 1 T (right panels) for m > 
(upper panels) and m < (lower panels) with R = 50 nm. 



In the limit of t' >> m + /3 we obtain s 2 _ = (to + /3) 2 [(to + 
0) 2 — \]/(t') 2 and thus the low energy solutions are given 
by 



1 



±-v/(to + /3) 2 ((to + /?) 2 -1). 



(22) 



For bilayer graphene, the wavefunctions are eigenfunc- 
tions of the following operator: 



J z = L z + hr z + HS Z , 



where now 



-28 2 [(m + P) 2 - (t') 2 /2}+8 4 



0, 



(20) 



1 (-1 

2 I 







9 - - ( Ul 
bz ~ 2 I cr. 



(23) 
(24) 



where s = e — tu is the energy shifted by the average 
potential. In the next section we report the results for 
the case of Ux = — U2 = Uj, where, the average potential 
u is zero. In the limit <5 — > 0, the quartic equation is 
reduced to a quadratic equation in s 2 and we obtain the 
real solutions 

4 = (m + /3) 2 + {t') 2 /2 



are 4x4 matrices. 

In bilayer graphene the components of the current den- 
sity are given by 



]x = v F 



(25)" 

The angular current can be calculated from the following 
relation. 



± v / (i') 4 /4+(™ + /3) 2 (l + 2 , (21) 

which results in four solutions for the energy. These are 
real when |to + /3| > 1. In the opposite case of |m+/5| < 1 
(or equivalently — 1 + /3 < m < 1—0) we have s 2 _ < and 
consequently the corresponding energies are imaginary. 



v F 



t (°* v m 











(26) 



Where, £((fS) is given by Eq. ©. We obtain for the 
angular current in the K-valley, 

JK = V F ((/>c(j) D + (j)* D (j)C ~ <i>*A^>B ~ $*B<t>A), (27) 



(i 




FIG. 7: (Color online) (a) Electron energy levels of a graphene 
single layer quantum ring as function of external magnetic 
field Bo for the same parameters as used in Fig. 5. Black 
curve shows the ground state energy (b) Ground state expec- 
tation value of L z /h as function of magnetic field for both K 
(black dashed curve) and K' valley (black dash-dotted curve). 
Expectation value of S z /h versus magnetic field is plotted in 
upper inset for K-valley and in lower inset for A"'-valley. Blue 
solid curve shows the expectation value < J z > which is the 
same for both valleyes. 



and the total angular current is given by j = jx + 3k' 
Where, the four spinor components are: 



MR) = h 



MR) 
MR) 
MR) 



m + f3 
e — tu\ ' 

(e - tu-y) 2 - (m + /3)(m + (3 - 1) 
i'(e-rui) ' 
(m + p)[{e - t Ui ) 2 - (to + /3)(m + /3 - 1)] 



t'(e - TUi)(e - tu 2 ) 



(28) 



Note that the radial current can be calculated through 



Jr = v F 



^ I o t(M ]f 



(29) 



where, j r = iv F (4>* A <f)B - <t>* B <t>A + 4>* c 4>d - <t>* D M) = 
for the case of an ideal ring. Using Eq. (|2T|) . the total 
current density becomes, 



T = ±l 

2v F (m + f3) 
e — tui 



^ [(e - riii) 2 - (to + P){m + $ - 1)] 
t' 2 (e — tu 1 )(e — tu 2 ) 



2 1 




FIG. 8: The angular current density in the (a) K- valley, (b) 
if'-valley and (c) the total current density of a monolayer 
graphene quantum ring as function of external magnetic field 
Bq for the ground state energy shown by the black curve in 
Fig. Era). 
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FIG. 9: (Color online) Lowest energy levels of a bilayer 
graphene quantum ring as function of ring radius R with (a) 
Bo = T and (b) B = 5 T for U b = fOO meV and to- 
tal angular quantum number — 10 < m < — 1 (red curves), 
1 < m < 10 (blue curves) and m — (green curves). The 
insets are an enlargement of the small energy and small R 
region. 
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FIG. 10: (Color online) Lowest energy levels of a graphene 
bilayer quantum ring as function of total angular momentum 
label m for B = -5, 0, 5 T with U b = 100 meV and R = 50 



FIG. 12: 
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(Color online) The same as Fig. 111! but for U b — 
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FIG. 11: (Color online) Electron and hole energy states of a 
graphene bilayer quantum ring as function of external mag- 
netic field Bo for Ub = 100 meV and R — 50 nm. The energy 
levels are shown for the quantum numbers — 10 < m < — 1 
(red curves), 1 < m < 10 (blue curves) and m = (green 
curves) . 



(30) 



B. Results 



The dependence of the spectrum on the ring radius, for 
Bq = T (upper panel) and Bo = 5 T (lower panel) is 
shown in Fig. [5J for a gate potential Ub = 100 meV, which 



for 73 = T opens up a gap in the energy spectrum. As 
compared to the single layer quantum ring results of Fig. 
[U we find two main differences: i) for R — > there are two 
states inside the gap, and ii) we have a second set of levels 
that for large R are displaced in energy by t. In the limit 
R — > the most important term in the dispersion relation 
is (to + /3) 2 [(to + fi) 2 — l] . For m = — 1, 0, 1 the behavior 
of the spectrum is different and the corresponding energy 
levels don't diverge when R — > 0. The same behavior was 
found for the single layer results, but only for to = 0, — 1. 
Previously, we found that for rings with finite widths the 
spectrum exhibits anti-crossing points which arise due to 
the overlap of gate-confined and magnetically-confined 
states. In the present model the carriers motion along 
the radial direction is neglected and consequently we have 
level crossings instead of anti-crossing points in the spec- 
trum. The dependence of the energy eigenstates on 
the angular momentum index to is displayed in Fig. [10] 
for U b = 100 meV, R = 50 nm, with B = -5 T (di- 
amonds), T (circles) and 5 T (triangles). Due to the 
finite bias in this case, the fourth-order character of the 
dispersion Eq. (fT^j) causes the curves to exhibit a Mex- 
ican hat shape. The energy minima for Bq = —5,0,5 
T are respectively given by m = —1,-10,-20. In Fig. 
[HI the energy levels are plotted as function of magnetic 
field, for a quantum ring with Ub = 100 meV, R = 50 nm, 
and for —10 < to < — 1 (red curves), 1 < to < 10 (blue 
curves), and to = (green curves). These results are 
very similar to those found for a finite width ring and 
exhibit two local minima that are separated by a saddle 
point. In the case of finite width quantum rings there are 
additional energy levels correspond with states that are 
partly localized outside the ring. Figs. rriT a.b) show the 
asymmetry between the electron and hole states, caused 
by the bias. It is seen that the electron and hole energies 
are related by Eh(m,B ) = — E e (— to, — B ), where the 
indices h(e) refer to holes (electrons). In the absence of 
bias, the electron-hole symmetry is restored, as shown in 
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FIG. 13: (Color online) Lowest energy levels of a bilayer 
graphene quantum ring as function of the gate potential Ub 
when Bo = T (left panels) and Bo — 1 T (right panels) for 
m > (upper panels) and m < (Lower panels) with R — 50 
nm. 

Fig. 1121 for a ring with R = 50 nm and the parabolic 
energy spectrum is recovered with zero energy gap. 

In Fig. [T2J the energy branches are plotted as func- 
tion of the bias, for both the zero field case (left panels) 
and for Bq = 1 T (right panels), with m > (upper 
panels) and m < (lower panels). Notice that the fig- 
ures are quantitatively similar to those found previously 
for a quantum ring made of a single layer of graphene 
where the gate potential Ub has a similar effect as the 
mass term A. The differences are that for Bq = T the 
degeneracies are now: i) E(Q) = E(l) = E(—l) and ii) 
E(m) = E(—m) for \m\ > 1. In the presence of the mag- 
netic field a gap is opened even for Ub — meV, which is 
more clearly illustrated in the inset of the right-bottom 
panel of Fig. [T3] Notice that here we found that for 
m = — 1 and m = — 2 no real energy solution is found for 
Ub below some critical value. 

Figs. [LlTb.c) shows the ground state expectation value 
of the angular momentum versus the magnetic field to- 
gether with the quantum number m (blue solid curve) 
which is an eigenvalue of the total momentum operator 
J z . Notice that the expectation value of J 2 , i.e. m, is dif- 
ferent in the K and K' valley which was not the case for 
monolayer graphene. The energy levels for the K (solid 
red curves) and K' (dashed blue curves) valleys are de- 
picted in Fig. HUa). Black curve (gray curve) shows the 
ground state energy for the K- valley (K'- valley). Notice 
that in the considered case we find that the difference 
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FIG. 14: (Color online) (a) Electron energy levels of a bilayer 
graphene quantum ring as function of external magnetic field 
Bo for a quantum ring of radius R — 50 nm and with Ub = 
100 meV for both the if-valley (solid curves) and the K'- 
valley (dashed curves). Black curve shows the ground state 
energy of the energy spectrum in if -valley whereas the gray 
curve the corresponding ground state energy of the K'-valley 
(b) Ground state expectation values of L z /h, S z /h, r z /h as 
function of magnetic field in the _R"-valley. Blue solid curve 
shows the expectation value of J z /h operator, (c) The same 
as (b) but, for if'-valley. 



between < J z >= mh and < L z > is about (0.7 — 0.8)h 
for both K and K' valleys. 

The ground state angular current of a bilayer graphene 
as function of magnetic field Bq in the if -valley jx , the 
if'-valley jx' and the total angular current j is shown 
respectively in Fig [THa). (b) and (c). In the case of a 
bilayer graphene quantum ring, the energy levels in the 
vicinity of the K and K 1 points are different because of 
the valley splitting and consequently the total angular 
current versus magnetic field is a more complicated saw- 
tooth function. Notice that the angular current for the 
K or K 1 valley is not zero at Bq = which is due to 
the valley polarization whereas the total current is zero 
at B = 0. 



IV. SUMMARY AND CONCLUSIONS 

In summary we considered the behavior of carriers in 
single and bilayer graphene quantum rings within a toy 
model. Our approach leads to analytic expressions for 
the energy spectrum. In our simple model we are not 
faced with the disadvantages of the nature of edge effects 
which appears in quantum rings created by cutting the 
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FIG. 15: The ground state angular current density in the 
(a) A'-valley, (b) i^'-valley and (c) the total current density 
of a bilayer graphene quantum ring as function of external 
magnetic field Bo with Ub = 100 meV and R — 50 nm. 



layer of graphene (or lithography defined quantum rings). 

We found an interesting new behavior in the presence 
of a perpendicular magnetic field, which has no analogue 
in semiconductor-based quantum rings. In single layer 
graphene quantum rings only for A > hvp/2R we found 
the opening of a gap in the energy spectrum between 
the electron and hole states. For both single layer 
and bilayer graphene quantum rings the eigenvalues 
are not invariant under a Bq — > —Bq transformation 
and in the case of bilayer the spectra for a fixed total 
angular momentum index m, their field dependence is 
not parabolic, but exhibit two minima separated by a 
saddle point. The persistent current exhibits oscillations 
as function of the magnetic field with period Qq/ttR 2 
which are the well-known Aharonov-Bohm oscillations. 
Because of the valley splitting in the energy spectrum 
of bilayer graphene the total current density versus 
magnetic field is a more complicated sawtooth function. 
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